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Analysis of future space energy needs has indicated a variety of 
missions requiring total energies exceeding those of past missions by 
two or three orders of magnitude (1 to 100 kW) . A major part of such 
a system will be the power converter and inverter subsystems which 
process this energy to meet user requirements • 

One of the major components within such a power subsystem is the 
magnetic device. This component is used for voltage and current 
transformation* input/output filter sections* parameter sensing* and 
semiconductor control. Recently, attention has been given to the 
different loss characteristics of magnetic devices observed under 
various excitation conditions. 

This dissertation examines the basic magnetic properties under 
various operating conditions encountered in the state-of-the-art 
DC-AC/DC converters. Using a novel core excitation circuit* the basic 
B-H and loss characteristics of various core materials may be observed 
as a function of circuit configuration* frequency of operation* input 
voltage* and pulse-width modulation conditions. From this empirical 
data* a mathematical loss cbaracceristics equation is developed to 
analytically predict the specific core loss of several magnetic mate- 
rials under various waveform excitation conditions. 

Previous magnetic core characteristics were observed using either 
a voltage source or a current source. For a core being excited in 
either of these modes, the voltage or current characteristics are ob- 
served as the independent variable. However, the excitation circuit 
used for this dissertation uses a high power wide bandwidth operational 
amplifier which senses the derivative of flux with time as the con- 
trolled parameter. Under these conditions flux is observed as the 
independent variable with voltage and current as the dependent vari- 
able. 
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Power measurements were obtained using a true RMS wide bandwidth 
wattmeter. The B-H loop characteristics were observed using an X-Y 
oscilloscope. This experimental excitation set-up allowed the mag* 
netic characteristics to be measured and observed under square » ramp, 
and sinusoidal flux conditions. 

Using this flux controlled excitation, magnetic core character* 
istics were developed for constant values of ras^ flux and frequency. 
These characteristics will show the circuit designer for the first 
time the functional relati rnships between induction level and specific 
core loss as a function of the two key DC-DC converter operating 
parameters of input voltage and duty cycle. 

Core excitations involving symmetric and asymmetric conditions 
for both sinusoidal and square wave waveforms were investigated ana- 
lytically using three specific core loss calculation methods. Experi- 
mental data was used to compare the accuracy of these methods. The 
predicted and measured specific core losses were within 10 percent of 
each other for all conditions analyzed. This deviation is within the 
error band of the instrumentation. Therefore, these methods will ana- 
lytically predict the specific core loss for any type of voltage exci- 
tation encountered in DC-AC/DC converter configurations using only 
one set of characteristic curves for the specific magnetic material, 
core configuration, and tape thickness. 
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CHAPTER I 


INTRODUCTION 

Magnetic materials offer the power electronics designer a wide 
range of magnetic characteristics (Ref. 1) from which he must select 
those materials applicable to his particular DC-AC/DC converter design. 
Like their semiconductor counterparts , magnetic devices display fami- 
lies of operating characteristics v;hich are functions of frequency, 
temperature, basic material composition, fabrication techniques, and 
driving source impedance (Ref. 2). The interaction between these var- 
ious functions may be explained by ferromagnetism, material properties, 
or domain wall movement (Ref. 3). However, for power electronics appli- 
cations it is desirable to provide a set of magnetic characteristics 

which are easily understood and accurately predict the performance of a 
particular magnetic material when applied in a specific DC-AC/DC con- 
verter configurriion. 

Trio dissertation examines the B-H and loss chars-^teristics of 
several magnetic raatarials used lii present state-of-the-art power proc- 
essing circuits. The effects of variables such as frequency, voltage 
level, voltage waveforms, and modulation upon the inherent material 
characteristics are explored using only externally measurable quanti- 
ties . 

Traditionally derived core loss data are presented as a function 
of flux density with frequency as a parameter. In general, hard exper- 
imental data is Taken with sinusoidal voltage excitation of 60 or 
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400 hertz. The family of characterist.^ c core loss curves are then gen- 
erated by using a mathematical formula such as that proposed by Legg 
(Ref, 4). However, this method is only valid for low values of flux 
density. Only recently has attention been given to the different loss 
characteristics of magnetic materials observed between sinusoidal and 
square-wave excitation conditions at 10,000 hertz (Ref. 5). In the 
past, only circuits having similar excitation conditions rou^d utilize 
these loss characteristics. Newly available state-of-the- rt instvu- 
mentation equipment, and mathematical techniques developed in this dis- 
sertation allow the power electronics designer to calculate the specific 
core loss for various nonsinusoidal, and asymmetrical square wave exci- 
tation conditions using only these loss characteristics. 

The square wave B-H loop and associated core loss characteristics 
may be observed in a unique manner by means of an excitation source in 
wnicii the flux rate of change (d<(>/dt) is observed as the independent 
variable, and ttie power loss is measured using a wide bandwidth (DC- 
300 K Hz) wattmeter. This excitation circuit, which is fully described 
in Chapter II, allows the control of an additional variable, d<()/dt, 
rather than just the driving frequency as is the case for all previc* 
sinusoidal core loss measurements. If the excitation is restricted tv 
that of a constant frequency square wave voltage having a variable duty- 
cycle, the d(J)/uc effects may then be determined for various constant 
values of independent of frequency. For ,.ae first time, the dis- 

tribution of the B-H loop area is shown to be dependent upon the posi- 
tive and negative d(|»/dt ramp rates. The wideband wattmeter allows 
core loss data to be accurately measured under these conditions. This 
procedure allows a true separation of the variables d>/dt, frequency, 
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and ITils data is then used to develop the functional relation- 
ships of Hysteresis loss with and eddy current loss with and 

d(()/dt. 

Sinusoidal core loss characteristics for a particular magnetic 
material and lamination thickness are generally shown with flux density 
as the independent variable, specific core loss as the dependent vari- 
able, and frequency as a pa>-air. er (Ref. 6). Data collected with the 
new excitation circuit is based upon values of constant d(|)/dt which 
corresponds with square wave excitation. The core loss characteristics 
are stiown witn frequency as the independent variable, specific core 
loss as the dependent variable, and d^/dt and flux density as param- 
eters. 

In Chapter III these loss characteristics, the B-H loop behavior, 
and the instantaneous power waveshapes observed for the various wave- 
forms are used to develop three methods which analytically predict the 
specific core loss for magnetic materials. A pulse width modulated 
square wave witn a 10 to 50 percent duty cycle is then used to compare 
the experiment. and analvtical accuracy of each method. 

The first method, which is similar to the sinusoidal analysis 
tecnniques already in common use, allows the specific core loss to be 
calculated directly fro? the characteristics knowing the maximum flux 
density (B^^) and the frequency of operation. To calculate the value 
of for the excitation pulse, the cross sectional area of the core, 

number of turns, excitation voltage level, and pulse width must be 
known. In the second method, the specific c>re loss Is calculated by 
using d(^/dt and ai* ‘'equivalent** excitation frequency. 

The value of d'J)/dt is directly proportior.al to the input voltage 
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and inversely proportional to the number of turns on the core. The 
"equivalent" frequency is inversely proportional to twice the pulse 
width. The third method is developed as an extension of the second 
method. Tiie excitation pulse is subdivided into a discrete number of 
time steps. The average voltage level is determined for each time 
step. This voltage level along with the number of turns on the core 
correspond to an incremental value of d0/dt. This :orresponds to an 
incremental loss at some "equivalent" frequency. The total specific 
core loss is then the summation of these incremental losse** divided 
by total period. These methods cover the functional relationship 
between ir&duction level and specific core loss as a function of exci- 
tation voltage level and duty-cycle characteristics. 

The concepts developed for the square wave in the third method 
are extended to analyze other types of v:aveforms in Chapter IV. A 
general equation form is used to relate the specific core loss to 
d(^/dt, and the equivalent frequency of the excitation. The coeffi- 
cient and expedients that relate these variables are dependent on the 
magnetic material, and the tape thickness, or grain size. Using this 
general equation a computer algorithm is then developed which auto- 
mates the analysis procedure. Analytical and experimental results 
are compared for sinusoidal and nonsinusoidal examples using three 
different core materials and configurations. As a final example the 
specific core loss characteristics for 1/2 mil Supermalloy mat^^rial 
under sinusoidal excitation are generated using the characteristic 
equation for square wave excitation and the computer algorithm. 

Finally, Chapter V presents the conclusions of this work and sug- 
gestions for future research in this area. Spec tic core loss data 
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for several magnetic tape materials and ferrite are presented in 
Appendix A. A glossary i>f svmbi>ls and their associated units are 
in Appendix B. 



CHATTER II 


INVESTIGATION OF MAGNETIC CORE CHARACTERISTICS 
2.1 Introduction 

The analytical core loss calculation methods developed in this 
dissertation are based upon the core characteristics observed under 
controlled flux excitation. The width of the B-H loop was found to 
be a function of both the flux rate of change (d(t>/dt) and the maximum 
flux density (Bj^) . 

In all the specific core loss (SCL) characteristics presented, 
both positive and negative values of dcfr/dt were held constant over 
one cycle. These d<()/dt conditions allow symmetrical square wave 
voltage excitation of the magnetic material under test. With the data 
recorded in this manner Faraday’s Law (Ref. 7) is then used to calcu- 
late the d(j)/dt value directly using the excitation voltage level, 
and the number of turns on the core. In an actual DC-DC converter 
application the excitation voltage level is the input voltage being 
applied to a particular magnetic device. Using the SCL characteristics 
as developed by this method, the power electronics designer may for the 
first time calculate the SCL for the magnetic device directly knowing 
or4ly the input voltage, duty cycle, and the frequency of operation. 

If the effective cross sectional area (A^) is known, then the 
flux density in the magnetic material may be calculated. To accurately 
determine the value of A^ the physical cross sectional area of the 
core (A^>) and the stacking factor (SF) must be known. Four methods of 


6 



7 


determining the SF value for cut cores are presented in this section. 
The results of these methods indicate that the SF value for this 
core configuration can be measured to within 5 percent of its actual 
value. 

Tlie excitation circuit allows three different source configura- 
tions to be simulated. By using voltage feedback the excitation source 
represents a voltage source. If current feedback is used, the excita- 
tion source represents a current source. When flux feedback is used, 
the observed magnetic characteristics become independent of source im- 
pedance. Thus, true sinusoidal or square wave voltage excitation can 
be easily generated. This excitation circuit eliminates many of the 
problems encountered with conventional square wave excitation sources. 

Tlie instrumentation system was designed to efficiently collect 
the ^ore loss data required for the development of the desired SCL 
characteristics. The system was also used to run the desirea experi- 
mental waveshapes against which the analytical core ss calculations 
were then compared. The instrumentation system utilizes a Phillips 
PM 3252 multiplying oscilloscop:* for observing instantaneous p(t) 
waveshapes. A Tektronix 7704 mail: frame is used to simultaneously ob- 
serve e(t), i(t), p(t), 4>(t), and the instantaneous B-H loop. 

2.2 Review of Magnetics Equations 

All available presentations of specific core loss characteristics 
display flux density as the independent variable, SCL as the dependent 
variable, and frequency as a parameter. To use this data requires that 
the A^ of the core be accurately known. This parameter is not always 
readily available and is inaccurate in most of the published data. 
However, for the majority of power electronic applications the input 
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voltage 9 frequency of operation, number of turns, and duty cycle con-* 
ditlons of a magnetic device are known. 

In a magnetic device a time-varying magnetic field produces a 
voltage which may establish a current in a suitable closed circuit. 
Faraday's law (Ref. 7), expressed in the C6S electromagnetic system 
of units, shows the relationship between this changing magnetic field 
and the induced voltage: 

e(' ) * -(d(^/dt)xlO~® volts (2.1) 

The minus sign indicates that the voltage is in such a direction as to 
produce a current whose flux, if added to the original flux, would re-- 
duce the magnitude of the voltage. This statement that the induced 
voltage acts to produce an opposing flux is known as "Lenz*s law." 

If the magnetic field is coupled by an N turn conductor, then 

e(t) “ -N(d(^/dt)xlO*® volts (2.2) 

By use of equation (2.2), transformer equations for both sinusoidal 
and quare wave voltage excitation may be derived. If e(0 represents 
a square wave of voltage, then ((>( ) would be a ramp function with its 
slope proportional to the volts/tum ratio. Over one half cycle 
(0 ^ t j< T/2), e(: ) has a constant value E. For the same interval 
^i‘c) gc *s from -<()^ to The standard form of the transformer 

p 'Li ition for square wave excitation is obtained by integrating equa- 
tion (2.2) over this half cycle interval: 
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T/2 

0 



(2.3) 


(2.4) 


f » t'^ Hz (2.5) 

Substituting (2.4) and (2.5) into (2.3); 

E = -4Bj^^SFf Nxio"® volts (2.6) 

If e(t) represents a sinusoidal voltage, then (p(t) would be a 
cosine function. Over one half cycle (0 <. t <. n/ai) , e(t) is a sine 
wave with an an^litude of E. For the same interval 4»(t) goes from 
to The standard form of the transformer equation for sinus- 

oidal excitation is obtained by integrating equation (2.2) over this 
Interval: 





ir/ui 

0 
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2<j> 


Ex 10® 


M 03N 


Rearranging terms 


E = (I)j^ Nxl0~® volts 


-8 


(2.7) 


Since 

(jj = 2irf rad-sec ^ (2.8) 

Substituting (2.4) and (2.8) into (2.7) 

E = SF 2ti£ NxlO”® volts 

or 

E * -4.44 SF f nxlO ^ volts (2.9) 

By using Equation (2.9) the flux density for a given excitation 
condition may be calculated. Once this value is known, published SCL 
data for sinusoidal excitation may be used to directly determine Che 
SCL for a particular magnetic material. To the author *s knowledge, 
published data for square wave excitation conditions is essentially 
nonexistent. However, using the methods described in this disserta- 
tion square wave characteristics may be easily obtained. The square- 
wave SCL characteristics have a distinct advantage over the sine SCL 
characteristics because the effects of an additional parameter d(()/dt 
may also be examined using the same set of characteristics. In addi- 
tion, by rearranging the representative axes with frequency as the in- 
dependent variable, and SCL as the dependent variable, with Bw and 

1 1 

d(J)/dt as parameters, the characteristics become easier to use. In 
Chapter 4 the square wave SCL characteristics will be used to gen- 
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erate the sinusoidal SCL characteristics directly. 

The SCL characteristics shown in Appendix A represent the total 
core loss as a function of frequency. If the energy loss is examined 
on a per cycle basis (Ref. 8), the effects of d(J>/dt, and frequency 

may be observed. The energy loss associated with the hysteresis loop 
is readily determined using Warburg's law (Ref. 9): 




dB 


watt-sec 


( 2 . 10 ) 


In the instrumentation system used, the input current to the core 

under test was observed as the abscissa, and the integral of the input 

voltage was observed as the ordinate. Using Equation (2.2), and the 

chain rule. Equation (2.10) may be rewritten in terms of the input 

voltage 6^(^) and the input current i^(t): 

-d:“T 


W, 


4tt 



"f 


t*0 


4it 



= -(4ttLA 
h e 


-T 

|nl“^ l^(t)JfA“V^ej^(t)jdt 

-1 

) / p^(t)dt watt-sec 

t=0 


( 2 . 11 ) 


Equation (2.11) indicates that the area enclosed by the observed 
B-H loop is proportional to the energy (hysteresis) loss per cycle. 

Separation of the eddy current, hysteresis, and anomalous loss 
components (Refs. 10 to 13) of the total core loss has been explored 
using only sinusoidal excitation methods. However, by using pulse 
width modulated square wave excitation techniques a direct separa- 
tion of loss variables may be accomplished. 
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2.3 Calculation of Stacking Factor 
Four different methods were used to determine the stacking factor 
of the cut cores used for this dissertation. The oil displacement 
method and t^e cross sectional scan method used the physical properties 
of the core to determine the stacking factor value. The sinusoidal 
and square wave excitation methods used the electrical properties of 
the core, and an assumed value for the maximum flux density. Super- 
malloy cut cores with 1/2, 1, 2, and 4 mil tapes were used and the 
stacking factors obtained using these four methods were compared with 
published values. Table I shows the results of this comparison. The 
four methods used are discussed in detail below: 

Displacement method . - In this method Archimedes* principle 
(Ref. 14) is used to determine the true d^ placement volume of the 
core. First, the core is weighed in air; 

W * p V- 

core . core Fe 

air 

or 

Vp = w (p )~^ 

Fe core . core 
air 

TABLE I. - COMPARISON OF FOUR METHODS FOR DETER^^INING 



STACKING 

FACTOR WITH 

PUBLISHED 

VALUES 


Tape 

thickness , 
mils 

Displacement 

method 

Sine ex. 
method 

Square ex. 
method 

Cross sec. 
method 

Published 

value 

0.5 

0.865 

0.885 

0.843 

0.914 

0.80 

1.0 

.888 

.928 

.856 

.895 

.90 

2.0 

.932 

.965 

.934 

.941 

.92 

4.0 

.927 

i.OO 

.934 

.967 

.94 
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Next the core is weig^ed in an oil of known density* The displace- 
ment weight AW is then 

AW»W -W =p.,V 

core . core . . oil core 
air oil 


or 


V = AW(p 
core oil 


-1 


The stacking factor is the ratio of the air volume of the core to the 
oil volume of the core: 


Stacking factor = 

Example: 

Core: 1/2 mil Supermalloy cut core 


i?e 


core 


p =8.7 gm-cm“^ 
core ® 


W. 


core 


. : 345 grams 


"air 

Wco reoil = ^^2.5 grams 


■* 0.928 gm-cm 
AW = 42.5 grams 


-3 


Stacking factor = (345/8.7)^(42.5/0.928) = 0.865 

Excitation methods . - These are indirect methods in which the 
saturation flux density is used to determine A^. The core is excited 
until the current waveform indicates that the core is in the satura- 
tion region. Then the excitation voltage level, and frequency are 
measured. Using the transformer equations derived in Section 2.2, 
may be calculated: 

For the sine wave: 

-8 

A = E(4.44 B,,f NxlO 
e ^ M 


For the square wave: 


Ag » E(4B^f NxlO"®) 


-1 
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The stacking factor is expressed as 


Stacking factor 


A 

Vc 


These values are measured at several frequencies and an average stack- 


ing factor value is calculated. 


Example : 



Core: 1/2 mil Supermalloy 

A : 
c 

2.277 cm^ 

i') : 10 turns 


7.2 KGauss 

For sine wave data: 



f. 

e, 

Ae» 

hertz 

volts 

cm^ 

600 

3.86 

2.01 

700 

4.48 

2.00 

800 

5.13 

2.00 

900 

5.83 

2.02 

1000 

6.5 

2.03 



A. = 2.01 cm‘ 
e 


Stacking factor ■ 2.01/2.277 « 0.883 


For square wave data: 

i, e, 


hertz 

volts 

cm^ 

600 

3.36 

1.94 

700 

3.86 

1.91 

800 

4.41 

1.91 

900 

5.07 

1.95 

1000 

5.44 

1.89 



• 1.92 cm' 


Stacking factor * 1 . 9212,277 ■ 0.843 


Cross sectional scan . - In this method a ml cropho cograph of the 
core cross section Is used to determine the stacking factor value. An 
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optical micrometer is used to measure Che adhesive thickness and the 
actual tape thickness. Since the tape width is constant » the stacking 
factor is the ratio of the tape thickness to Che tape plus adhesive 
thickness width or 

Stacking factor ■ T /(T + T „ , ) 

® tape tape adhesive 

If these values are measured at several points across the core sur- 
face, and average stacking factor value may be calculated. In Fig- 
ure 1 the microphotographs for the four tape thicknesses are shown. 
Example : 

Core: 1/2 mil Supermalloy 


Tape thickness, 
in. 

0.000609 

.000578 

.000537 


Adhesive thickness, 
in. 

0.00003125 

.0000375 

.0000375 


0.000546 0.0000468 

= 0.000546 (in.) T , - 0.00005125 (in.) 

tape adhesive 

Stacking factor = 0.000546/(0.000546 + 0.00005125) 

= 0.914 

2.4 General Experimental Set-Up 

Previous excitation methods have used feedback control schemes 
(Ref. 15) to reduce the harmonic content of the sinusoidal excitation 
by controlling the apparent source impedance. For this work, feedback 
control is used to achieve the desired voltage excitation. If a square 
wave voltage is required, a triangular reference voltage is used. If 
an impulse voltage is required, a square wave reference voltage is used. 



(an/2mil Supermalloytape. 



(b)l mi! Supermallf'v tape. 

Figure i.- Cross section Microphotogr?'‘* of 1/2,1, 
and 4 mil Supermalloy cut cores. 
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Aad» if sinusoidal voltage excitation is required, a cosinusoidal ref- 
erence is used* Figure 2 shows the voltage, current, and flux relation- 
snips for each of these flux reference conditions . 

£he complete instrumentation system is shown in Fi^ ire 3. A high 
power operational anq>lifier was designed which allowed option of 
selecting eitner voltage, current, or the integral of voltage (flux) 
as the feedback, parameter. A Wavetek model 146 signal generator was 
used to generate the desired symmetrical flux waveshape. An ICL 8038 
multifunction integrated circuit was used to generate the asymmetrical 
square and sinusoidal reference waveforms. The excitation frequency 
was measured with a Hewlett-Packard 5212 A electronic counter. A 
Clarke-Hess model 223 wideband wattmeter (0 to 300 kHz) was used tc 
measure the average input power into the core under test. This in- 
strument allows the power delivered by a 100 kHz square wave to be 
normally ir^.asured to within 2 percent. The internal current shunt 
in the wattmeter displays the input current signal used for several 
measurements. A Phillips PM 3232 multiplying oscilloscope was also 
used to obtain the instantaneous p^(t) in the core. The input volt- 
age to the core was measured with a Hewlett-Packard 3403C true RMS 
voltmeter. The Tektronix 7704 oscilloscope was used to observe i^(t), 
©iCt), (^(t), P^(t) as well as the instantaneous B-H loop of the core 
under test. 

2.3 Source Impedance Effects 

Source inpedauce effects discussed in the literature (Refs. 16 
and 17) show that the width of the B-H loop to be dependent upon the 
driving source impedance. By changing the feedback parameter on the 
operational amplifier, a voltagt. ^ R^^) or a current source 



e(t> -0.5Vmiv 


itt) -asA/oiv 


(pm »1050MaxfDiv 
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(R >> R, ) may be simulated. Tlie flux, current, and B-H loop charac- 
S la 

teristics for these two source configurations are suown in Figures 4 
and 5. For the same conditions, the width of the B-H loop is 

wider if the core is being driven by a current source. However, if 
the flux is observed for these two source conditions it can be seen 
that for sinusoidal current conditions, the flux approaches a square 
wave. This is a condition having a higher d4»/dt rate in the core 
than that realized under sinusoidal voltage conditions. Therefore, 
the observed B-H width is really a function of d({)/dt rather than 
being dependent upon source impedance. 









CHAPTER III 


ANALYTICAL CORE LOSS CALCULATICWS 
3.1 Introduction 

Analytical prediction of core loss in magnetic materials is 
important from both an economic (Ref« 18) and an engineering design 
(Ref. 19) standpoint. Accurate prediction allows design tradeoffs to 
be made between the areas of weight, efficiency, and performance. In 
addition new circuit developments employing magnetic devices may be 
more fully analyzed without hardware fabrication. Early work (Ref. 20) 
at frequencies of up to 200 kHz was aimed at measuring the magnetic 
properties of Iron. Separation of core loss into its specific compo- 
nents of eddy current, static hysteresis, and anomalous loss has been 
examined in great depth (Refs. 21 to 23) for low frequency (<1 kHz) 
sinusoidal excitations. 

With the growth of power electronics to higher power levels and 
the employment of new converter techniques, the need to understand the 
properties of magnetic materials under nonsinusoidal conditions is 
arising. Recent comparisons between sinusoidal and square wave core 
loss (Ref. 3) indicate a lack of theory to predict the loss difference 
observed between these two excitations. However, if the dependence of 
the characteristics upon the flux conditions described in Section 
2.3 IS expanded, the loss characteristics for these two excitations may 
be qualitatively observed, and predicted using easily obtained data. 

In this chapter a comparison between sinusoidal and square wave 
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covs loB8 is made, lha results are shown to be consistent with ^he 
values published in the literature. Next, pulse width andulatec s<;uare 
wave excitation is used as an analytical tool to achieve a true separa- 
tion between the d^/dt and effects. Finally, three oetitods ate 
developed to analytically predict the SCL for square wave excitation. 
These results are cospared with experimentally measured results ob- 
tained under the same ewitation conditions. 

3.2 Specific Core Loss Measurements 
To generate the SCL characteristics under square wave excltatiim 
for the magnetic materials studied, data was tdcen by holding d$/dt 
constant and varying the excitation frequency. The power input to the 
core was measured at regular frequency intervals. The actual values 
for SCL in wa'ts/pound and B{| in gauss were calculated from the 
raw data using the following formulas: 


Measured core loss in watts __ 

T7 : T~r — — B SCL(W/lb) 

Measurec core wei^t in pounds 


(3.1) 


Measured e in volts 

B„(gauss) (3.2) 

kA SFf NxlO"® 
c 


fdiere 


k = 4 for square wave excitation 


<= 4.44 for sinusoidal excitation 


Figure 6 shows the specific core loss for 1/2 mil cut Supermalloy 
under both sinusoidal (Fig. 6(a)) and square wave (Fig. 6(b)). In 
addition, comparisons are made with published sinusoidal data for these 
two excitation conditions. The major difference between the published 
and measured values for sinusoidal excitation may be accounted for by 
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several factors: 

Probable accuracy of the published data 
Accuracy of the core loss maasurement 
Limited raw core loss data available 
Driving source impedance 

Measurement of actiial excitation voltage level 
Value assumed for stacking factor 
Batch variations for the magnetic material 
Actual lamination thickness 

It is interesting to note that for equal values of the SCL 

for sinusoidal excitation is greater than for square wave excitation. 
Example: 

B^ a 2k gauss frequency ° 1 kHz 

SCL (sine) <» 0.078 W/lb 
SCL (square) = 0.058 W/lb 

From these results it can be seen that the SCL for square wave 
excitation is approximately 26 percent less than the SCL for sinus- 
oidal excitation, 'fliis result is consistent with previous observations 
(Ref. S) which compared loss conditions for these two excitations. 

To aid in the de>ralopment of methods for specific loss calcula- 
tions , data is replotted with the frequency and flux density axes 
Interchanged. Figure 7(a) shows SCL versus frequency with B^ as 
a parameter for cut 1/2 mil Supermalloy under sinusoidal excitation. 
Figure 7(b) shows these characteristics for square wave excitation. 

For square wave excitation conditions, d^/dt was held constant as well 
as frequency. This allowed the additional parameter of d^/dt to be 



Specific core loss, waWpound 



(a) Specific core loss characteristics for cut 1/? mil Super- 6)) Specific core loss characteristics for cut 1/? mil Super- 
malloy for sinusoidal excitation. malloy for square wave citation with B|y| and dfnfdt 

as parameters. 


Figure 7. - SCL characteristics for sinusoidal and square wave excitation with axes interchanged. 
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plotted on the same SCL characteristic (Fig. 7(b)). 

3.3 Separation of Loss Variables at Constant Frequency 
If excitation is restricted to that of a constant freauency pulse 
width modulated square wave, the d^/dt effects may be determined for 
various constant values of B^. Figure 8 shows the Instantaneous wave- 
forms for p(t), l(t), and e(t). Positive l(t) values correspond to 
H (Nl(t)L~^) values to the left of the B-axis while negative i(t) val- 
ues correspond to H values on the right of the B-axls. Also, p(t) 
values are both positive and negative which corresponds to energy sup- 
plied to or energy returned from the core. 

Figure 9 shows the instantaneous B-U loop of 1/2 mil Supermalloy 
measured for various d(|)/dt (change in duty cycle) values while holding 
c>'nstant at 6.76 K gauss, and frequency constant at 2 kHz. From 
this figure, it may be seen that as d4>/dt increases the B-H loop 
area for that portion of the cycle also Increases. Since the B-H 
area is proportional to t a energy loss, and the additional B-H area 
is proportional to d^/dt, it is assumed that the additional half cycle 
energy loss observed is proportional to d(|)/dt. 

This procedure was then repeated for several values of B^ with 
frequency held constant with the results shown in Figure 10. From this 
figure it is seen that for constant d(|>/dt, the total core loss in- 
creases with B^. If the lines of constant B^ could be linearly 
extrapolated, the ordinate (d(|)/dt « 0) Intercept would indicate the 
nysteresis increase with B^. The eddy current loss is dependent upon 
both dift/dt and B^. However, it is interesting to note that for a 
constant value of B^, the half cycle loss is linearly dependent on 
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1 A-T/Div 



0.2WDIV 


(a) B-H loop and Instantaneous p(t). 


• somv . mm 

/. m V 


p(t) 

0.2W/DIV 


irmm 7 ^" 


attf 

2 V/Div 


lb) Instantaneous pit), lit), and ert). 

Figure 8. -• Comparison of the four loss regions for the 
instantaneous B-H loop, pit), Mt), andeith 
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d<^/dt* This procedure shows that a true separation of the variables 
d$/dt, and may be achieved using pulse width modulated square 

wave excitation techniques. These characteristics will be used to 
qualitatively explain the differences observed in the B-H loop tor 
sinusoidal and square wave excitation. 

A comparison of the B-H loop characteristics under both constant 
(square wave (fig. 11(a)) and variable (sinusoidal (fig. 11(b)) d(j)/dt 
conditions measured while holding both frequency and B^ (relates to 
hysteresis loss) constant is shown in Figure 11. Figure 11(c) illus- 
trates some interesting conditions which, to the author*s knowledge, 
have never been shown in the literature. It has been shown (Fig. 10) 
that the total half cycle core loss is composed of a hysteresis loss 
which is proportional to and an induced eddy current loss propor- 

tional to both d(f)/dt and B^. 

Since B^^ is held constant, hysteresis loss is constant, and eddy 
current loss is dependent on d4>/dt only. Tlie following qualitative 
observations may be made. At point 1, d({»/dt for the sine wave is 
equal to zero while d<J>/dt for the square wave is equal to some posi- 
tive constant value greater than zero. Therefore, the core loss for 
the square wave is greater (wider B-H width) than that for the sine 
wave at point 1. As the flux density increases to point 2, d^^/dt for 
the sine wave reaches its maximum value. This value of d^i/dt is 
'greater than the d4>/dt for the square wave, and the corresponding core 
loss (B-H width) is greater. As the flux density increases to point 3, 
d(j)/dt*s for the sine and square wave are approximately equal. This cor- 
responds to an equal B-H width. As the flux approaches point 4 the 
same conditions as point 1 are repeated. 
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3,4 Analytical Calculation Methods 
Using the SCL characteristics for square wave excitation three 
methods were developed which analytically predict the SCL of a 
netic material tor pulse width modulated square wave excitation* liiese , 
Bietuods are discussed in detail below. 

Method i - Specific core loss based on of the core. '- The 

SCL characteristics presented in Figure ?(b) relate the total specific 
loss (hysteresis ^ eddy current » total loss) to frequency with B%| as 
a parameter. This data will be used to determine the core loss charac- 
teriscics witii pulse width modulated square wave core excitation. Ihis 
niethod is the most indirect approach to calculating the specific core 
loss. It requires an accurate measure of the effective cross sectional 
area of tt 'ire. Generally* this is a par aii^ ter whlciv;:^is: net .aceu*^ 
rateiy k. the power electronics designer. 

Figure r irhows the oscilloscope traces of a 1000 Hz square wave 
operating over d duty cycle range of iO to 30 percent. For each duty 
cycle the voltage for the positive cycle is adjusted to a maximum value 
of 4 volts. The DC offset in the core is held to a minimum.. Tlierefore, 
cue negative and positive half cycle volt-seconds are equal. 

In order to use the Figure 7(b) characteristics the value for 

each duty cycle condition must be determined. This value may be ob- 


tained as follows: 
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3.4 Analytical Calculation ^fethods 
Using the SCL characteristics for square wave excitation three 
methods were developed which analytically predict the SCL of a mag- 
netic material tor pulse width modulated square wave excitation. These 
methods are discussed in detail below. 

Method 1 - Specific core loss based on of the core. - The 

SCL characteristics presented in Figure 7(b) relate the total specific 
loss (hysteresis + eddy current * total loss) to frequency with as 

a parameter. This data will be used to determine the core loss charac- 
teristics with pulse width modulated square wave core excitation. This 
metiiod is the most indirect approach to calculating the specific core 
loss. It requires an accurate measure of the effective cross sectional 
area of ti “^re. Generally, this is a parameter which is net acci’- 
rately K.* the power electronics designer. 

Figure i thows the oscilloscope traces of a 1090 Hz square wave 
operating over a duty cycle range of 10 to 50 percent. For each duty 
cycle the voltage for the positive cycle is adjusted to a maximum value 
of 4 volts. The DC offset in the core is held to a minimum. Therefore, 
the negative and positive half cycle volt-seconds are equal. 

In order to use the Figure 7(b) characteristics the value for 

each duty cycle condition must be determined. This value may be ob- 


tained as follows : 
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From Faraday's law 

e - -N d(j)/dtxl0"8 (volts) (3.3) 

c 

Tlierefore 

—1 8 

d(t)/dt » N xio (maxwells /second) (3.4) 

For each cycle portion 

Ad» * (d(t/dt)At (3.5) 

AtxlO® (3.6) 

and, 

= A*(2Ag)“^ - -0.5 Ac(NA^)"^xlO® (gauss) (3.7) 

where: 

2 

« A X SF (cm ) 
e c 

N « number of turns 

e^, At related to duty cycle 

Using 1/2 mil Supermalloy as an example, for a frequency of 
1000 Hz, the values corresponding to each duty cycle width can be 

determined as follows : 

Bjj (positive cycle) - -0.5 (3.9) (100xl0“^) (10xl.97)"^xlO® 

“ -990 gauss 

(negative cycle) = -0.5 (-0.44) (900xl0~^) (10xl.97)"^xl0® 

* 1000 gauss 

(positive cycle) ; (negative cycle) 

Considering the positive half cycle only this calculation can be 
repeated for the remaining four cases with the results shown below: 
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Duty cycle, 
percent 

At-ys 

+ 

®c» 

volts 

4 * 

gauss 

10 

100 

3.9 

990 

20 

200 

3.9 

1980 

30 

300 

3.9 

2470 

40 

400 

3.8 

3858 

50 

500 

3.4 

4315 


This data can now be used with Figure 7(b) to determine the specif ic 
core loss associated with each duty cycle condition. The results of 
these calculations are shown on Figure 13. 

During the experiment the specific core loss was measured for 
each duty-cycle condition. A comparison between the measured and cal- 
culated results is shown in Figure 14* At the 50 percent duty-cycle 
condition the error between the measured and calculated specific core 
loss was: 

» 0.31 W/lb 

SCL = 0.339 W/lb 
meas 


- SCL 


percent error = 


X 100 


(3.8) 


0 339 - 0 31 

percent error » * ^ 339 — ^ percent 

From these results it may be seen that the specific core loss can 
be easily calculated for pulse width modulated square waves if the base 
frequency » excitation voltage level, duty-cycle, number of turns, and 
the effective cross sectional area of the core is known. Also, for a 
constant input voltage as the duty cycle increases, increases, and 
the specific core loss of the magnetic device increases. 
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Figure 13. - Specific core ioss characteristics for cut 1/2 roil Superroaitoy under duty 
cycle conditions. 
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Figure 14. - Comparison of caicuiated vs measured vaiues of SCL for 
1/2 mii Supermalloy as a function of duty cycie for a 1 kHz - 3. 9 
voit puise voitage. 
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Method 2 - Specific core loss based on d<^/dt In the core ^ - The 
specific core loss characteristics presented in Figure 7(b) relate the 
total specific core loss to frequency with and d^/dt as param- 
eters. If the parameter is ignored on this figure, the specific 

core loss can be calculated by using the d<^/dt parameter and an 
equivalent excitation frequency. 

In the discussion of Section 3.3 the core loss was shown to be a 
function of B^ and d<^/dt. For a pulse width modulated square wave 
with a constant positive d^/dt and constant negative d<b/dt, is 

n 

determined directly by the length of time that the pulse is on. The 
width of this "on pulse" has an "equivalent" frequency of 

f = = (2t )"^ (3.9) 

eq eq ' on 

Also, it was shown that the total specific core loss per cycle 
could be broken down into a half cycle loss dependent on d^/dt. 
Therefore, the total specific core loss may be determined using the 
following relationship. 


Total SCL 


(SCL/cycle) (excitation frequency) 


SCL 


eq 


1 SCL~ \ 

^ f‘ / 
eq / 


ex 


where 

SCL^ positive cycle specific core loss 

SCL negative cycle specific core loss 

f^^ positive cycle equivalent frequency 

f^q negative cycle equivalent frequency 

f 

ex 


(3.10) 


excitation frequency 
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Example : 

10 percent duty-cycle case: 

= 3.9 V 
c 

e” “ 0.44 V 
c 

£ ^ «= 1000 Hz 

C A 

t^ «= lOOxlO"^ sec 
on 

t“ ® 900x10”^ sec 
on 

Then 

f'*'- « (2t'^ = (2x100x10"^)"^ - 5 kHz 

eq on 

- - “1 ft -1 

f = (2t ) = (2x900x10"®) * 555 Hz 

eq on 

using Faraday’s la» 

d^/dt = e^ N~^xiO^ 

(d<(./dt)'‘‘ - 3.9(10)”^xio® - 39x10^ maxwells/sec 

(d^/dt)~ ■ 0.44(10)”^xl0® a 4.4x10^ maxwells/sec 

Replotting Figure 7(b) as Figure 15 with (d$/dt)^, (d(()/dt)”, 
and f” shown, the corresponding SCL’s may be read directly from the 
figure : 

SCL' = 0.082 watt /pound 
SCL” “ 0.0063 watt/pound 

Using the equation for the total SCL the following results are 
obtained: 

ro»isa.i(|M.2^)zooo 


Total SCL « 0.014 watt/pound 



43 



44 


summarizing the results for the other conditions we have: 


f 


base 


= 1000 Hz 


Duty cycle, 
percent 

eq 

Hz 

f 

eq, 

Hz 

dt 

6 

xlO max/sec 

^10^ max/sec 
dt 

10 

5000 

555 


39 

4.4 

20 

2500 

625 


39 

10 

30 

1670 

714 


39 

16 

40 

1250 

833 


38 

26 

50 

1000 

1000 


34 

34 

From Figure 

15 the 

SCL'*' 

and 

SCL can be 

obtained. Again 


using the total SCL equation: 

20 percent duty cycle: 

Total SCL - j “ 0.05A W/lb 

30 percent duty cycle: 

Total SCL = j + ^7^) 1000 “ 0.124 W/lb 

40 percent duty cycle: 

Total SCL = 1^7^+ 

30 percent duty cycle: 

Total SCL - + 2^) 1000 = 0.3 W/lb 

In Figure 16 these calculated SCL's are comparea with the meas- 
ured SCL values for these duty-cycle conditions. The percent error 
is the same as obtained with method 1 (8.6 percent). 

Method 3 - Specific core loss based on incremental d(t/dt in 
the core # - The technique developed in method 2 allows determination 
of Lhr specific core loss based on the d(()/dt value for the entire 
positive or the entire negative cycle of the square wave. In Sec- 
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Rgure 16. - Comparison of calculated vs measured values of SCL as 
a function of duty-cycle for a 1 kHz - 3. 9 volt pulse voltage. 
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tion 3*3 the dependency of the B«U width on d0/dt was developed* 

For a particular frequency and level a family of curves could be 

developed for various d(^/dt values as shown in Figure 17. Using 
actual square waves observed in this experiment (Fig. 12) , it can be 
seen that for any interval the actual exciting voltage is not a con- 
st' value. This would indicate that the corresponding d^/dt values 
over the cycle are not of constant value. If the waveform is broken 
down into discrete sample intervals the corresponding locus of points 
on the B-H loop may be plotted as At is incremented. The corre- 
sponding B-U width would be a direct function of the pulse height 
for that interval. Figure 17(b) shows how the pulse height, d^/dt, and 
At would translate into a locus of points on the B-H loop. This 
approach may be used to explain the B-H loop shapes observed under 
various excitation conditions. 

Relating these various d4»/dt rates to loss at the positive 
equivalent frequency, the SCL characteristics of Figure 17(c) may 
be used to determine the incremental core loss. 

The loss per interval is calculated as follows: 

= n At? » (3.11) 

^ 2f 


where n = number of increments 


At;!’ = (3.12) 

^ 2nf 

The energy dissipated in At^ would be: 

+ 

SCL4 

AE’T = SCL"^ At| ^ (3.13) 

^ ^ ^ 2nf 

If At^ were made small enough, the difference between the incre- 
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(a) family of curves for various 
d^^dt values. 




(c)SCLcharacteristics for various dt9/dt 
values. 


Hgure 17. * Family of B-H loops and SCL characteristics based on as a 
parameter. 
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mental excitation and the actual excitation would approach zero a 
limit. 


The total energy over the positive cycle is simply the summation 
of the incremental &e'*''s: 

N 


+ ' + 1 ^ 
eT • AE - 


T . , 1 

i-1 


sa] 


W-sec 


2nf^ ^ -'-'i lb-cycle 


Similarly the total energy for the negative cycle is 


M 


m 




SCL 


W-sec 


2mf 


- L_J i lb-cycle 


i=l i-1 

and the average power would be: 


4+ e; 

SCL^ = specific power loss * 


W/lb 


(3.14) 


(3.15) 


(3.16) 


T + T 

Applying this technique to the 10 percent PUM square wave case we 
have: 


Positive cycle 


Interval 
(Atj^ * 20 ps) 

®c’ 

volts 

M . 

dt ’ 

10^ max/sec 

scl[, 

W/lb 

1 

3.4 


34 

0.063 

2 

4.0 


40 

.092 

3 

3.9 


39 

.087 

4 

3.9 


39 

.087 

5 

4.0 


40 

.092 


E 0.421 
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Negative cycle 


laterval 
(At^ » 100 ys) 

volts 

, 

dt ■ 

10^ max/sec 

SCL". 

W/lb 

1 

0.2 


2 

0.00092 

2 

.4 


4 

.0044 

3 

.4 


4 

.0044 

4 

.4 


4 

.0044 

5 

.44 


4.4 

.0055 

6 

.5 


5 

.0074 

7 

.5 


5 

.0074 

8 

.55 


5.5 

.0091 

9 

.55 


5.5 

.0091 


E 0.05262 


-T 2. X 5 X 5000 


^ r (0.421) = 8.42x10“^ W-sec 


lb -cycle 


eI = > (0.05262) = 5.267x10"^ 

T 2 X 9 X 55 j ' ^ lb-cycle 


SCI^ 


8.42x10*^ + 5.267xl0~^ 
100x10"^ + 900x10"^ 


= 0.014 W/lb 


For Che remaiaiag cases we have 
20 percent: 


10 

r~\ + 

; SCL. = 1.764 

Lj ^ 

i=l 

8 

^^SClT = 0.298 

Lj ^ 

i=l 


At^ = 20 us 


At. = 100 us 

X 


1.764 


0.298 


. 2 >10 >2500 2 . 8 x62ii . 

200x10 + eooxio"'’ 
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30 percent: 


15 


= 4.554 = 20 us 


i«l 

7 


^sa“ = 0. 

i=l 


621 At^ = 100 us 


4.554 


0.621 


SCI. = ^ -.A^9 2 .x 7 X 714 ^ ,53 

^ 300x10 + 700x10 


40 percent: 


■20 

V^SClT “ 8.204 Att * 20 us 
/ / 1 i 

i=l 


y^scL? = 1.] 


_.117 At, * 100 us 

i=l 


7.73 


1.117 


5 . 2 X 20 X 1250 ^ 2 x _ 6 x 833 ^ ^.266 W/lb 

400xl0‘® + 600x10“° 


50 percent: 


25 

= 9.03 AtT = 20 us 

L-/ ^ i 

i=l 


5 

^SCL- - 1. 
1=1 


66 At, = 100 us 


9.03 


1.66 


SCUj. 


2 X 25 X 1000 • 2 X 5 X IQOQ 
500x10"^ + 500x10“^ 


0.347 W/lb 


A comparison of this value with the measured results is shown in 
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Figure 18. The error using this technique is less than that of methods 
1 and 2. 

The instantaneous p(t) waveforms are shown in fxgure 19. A 
Phillips PM3252 multiplying oscilloscope was used to obtain p(t). in 
the top photo a resistive load was used to verify the accuracy of the 
multiplier. For the square wave the growth of p(t) is fairly linear 
with time, which supports the preliminary assumptions used in devel- 
oping this method. 
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Figure 18. - Comparison of calculated vs measured values of SCL as a 
function of duty cycle for a 1 kHz - 3.9 volt pulse voltage. 



Resistive load 


Uncut core 


Cut core 


iCSi; 



CHAPTER IV 


GENERAL CORE LOSS Ai>IALYSlS 
4.1 Introduction 

Analysis techniques such as harmonic analysis and digital simu- 
lation of eddy current losses have previously been used for prediction 
of core losses of distorted flux waveforms (Ref. 24). Eddy current 
correction factors (Ref. 25) have also been employed in estimating the 
core loss of sinusoidal excitation considering only the fundamental 
and one harmonic component. In this chapter method 3 of the previous 
chapter is used to analytically predict the specific core loss for sev- 
eral distorted waveform excitations. 

Since this method is based upon the summation of a large number of 
sample points, tne digital computer was used for high speed core loss 
analysis. In order to minimize the number of subroutines of the com- 
puter algorithm, a general characteristic equation was derived using 
experimental data. This equation relates SCL to the equivalent fre- 
qverit,' and d^/it. An alternate equation which relates SCL to equiv- 
■it'n*: frequency and was also derived. For analysis of the SCL 

j.i/rscteristic3 a specific magnetic material having a known tape or 
grain thicknesses, only a constant and two exponent values must be 
supplied in the initial data format. 

Using the square wave SCL characteristics and the computer algo- 
ritnm, analyses were performed for sinusoidal and nonslnusoidal exci- 
tation conditions using Supermalloy, Square Permalloy 80, and ferrite 
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materials* As a final illustration of the method, the sinusoidal core 
loss cnaracteristics of a 1/2 mil Supermalloy material operating at 
1000 Hs is measured experimentally and calculated analytically. 

4.2 Development of Computer Algorithms 
The SCL characteristics for various types of magnetic materials 
is presented in appendix A. Due to the logarithmic nature of the char- 
acteristics, the SCL equation for either the d(f)/dt parameter or the 
parameter would have the following form: 

SCL = (4.1) 

The exponent x could be either positive (let x » d) or negative 
(let X = u) depending upon the parameter being examined. 'Hie constant 
a is also a logarithmic function of either d(|)/dt or The equa- 

tion for the constant is as follows: 

a * g(d<)/dt)^ (for the d(^/dt parameter) (4.2) 

or 

a * parameter) (4.3) 

Combining Equation (4*1) with the appropriate constant equation, the 
following two equation forms may be used to determine the SCL for a 
specific frequency and or d(J>/dt level: 

For d4>/dt: 


SCL = g(d(()/dt)‘^ (4.4) 

eq 

For 

SCL = s(Bj^)^ (4.5) 

Values for the constants g and s, and for the exponents c, d. 


z, and )i for various magnetic materials were determined from the 
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experimental data shown In appendix A. These calculated values are 
shown in tables II and III. In Equation (4.5) at low flux densities, 
the d exponent was determined by Stelnmetz to be 1.6 for electrical 
sheet steel (Ref. 26). For higher flux densities the value of the 
exponent may be as high as 2.5. The values of the exponent as meas- 
ured In these experiments falls within this upper limit. The value ot 
tne exponent was found to be dependent upon both the magnetic material 
and the tape thickness. 


TABLE II. - TABULATION OF THE CONSTANT AND EXPONENT 
VALUES USED IN EQUATION (4.4) 

Table II. 1. - Supermalloy-Cut Core 


Constant 

Tape thickness, mils 

0.5 

1.0 

2.0 

4.0 

c 

2.270 

1.563 

1.981 

1.959 

g 

2.537x10"^^ 

3.819x10"^® 

3.240x10"^^ 

4.420x10"^^ 

u 

-1.000 

-0.811 

-0.500 

-0.345 


Table II. 2. - Supermalloy - Uncut 


Constant 

Tape thickness, mils 

0.5 

1.0 

2.0 

4.0 

c 

1.951 

2.000 

2.090 

2.040 

g 

7.48x10”^^ 

4.430x10"^^ 

3.259x10"^^ 

1.124x10"^^ 

u 

1 

-0.408 

i 

-0.825 

-0.561 

-0.297 
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T ble II. 3. - Square Permalloy 80 - Uncut 


Constant 

Tape thickness, mils 

0.5 

1.0 

2.0 

4.0 

c 

0.985 

1.713 

1.704 

1.739 

g 

1.845x10"^ 

2.759x10"^^ 

1.201x10“^^ 

9.202x10"^^ 

u 

-0.349 

-0.225 

-0.226 

-0 • 258 


Table II. 4. - Ferrite 


1 Constant 

Material 

3B7 

c 

2.64 

g 

1.544x10“^^ 

u 

-1.6 


TABLE III. - TABULATION OF THE CONSTANT AND EXPONENT 
VALUES USED IN EQUATION (4.5) 

Table III.l. - Supermalloy - Cut Core 


Constant 

Tape thickness » mils 

0.5 

1.0 

2.0 

4.0 

d 

1.210 

1.307 

1.467 

1.670 

z 

2.038 

2.084 

1.948 

1.948 

s 

2.876x10"^^ 

1.065x10"^^ 

7.99x10"^^ 

3.231x10"^^ 


Table III. 2. - Supermalloy - Uncut 


Constant 

Tape thickness , mils 

0.5 

1.0 

2.0 

4.0 

d 

1.230 

1.300 

1.53 

1.750 

z 

1.940 

2.000 

2.1i0 

2.020 

s 

4.297x10"^^ 

8.676x10"^^ 

9.578x10"^^ 

1.056x10’^^ 
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Table III. 3. - Square Permalloy 80 - Uncut 


Constant 

Tape thickness , mils 


0.5 

1.0 

2.0 

4.0 

d 

0.624 

1.456 

1.459 

1.613 

z 

.976 

1.696 

1.722 

1.884 

s 

1.718x10"^ 

7.054x10"^^ 

7.330x10“^2 

9.447x10“^^ 


Table III. 4. - Ferrite 


Constant 

Material 

3B7 

d 

1.04 

z 

2.67 

s 

3.024x10"^^ 


The use of Equations (4.4) and (4.5) Is not restricted to a com- 
puter algorithm. P-y using a simple hand calculator and the analysis 
methods descrlb d In Section 3.4 the SCL for a particular magnetic 
material and tape thickness may be easily determined. An example of 
this process Is shovm In Figure 20. For the waveform shown the posi- 
tive half cycle frequency Is 8.77 kHz, and the negative half cycle 
frequency Is 12.5 kHz. The measured SCL for this waveshape Is 0.144 
watts /pound. If the SCL Is calculated using 5 ps Intervals, the 
calculated value Is 0.132 watt/pound. This represents an error of 
8.3 percent. This error will be reduced by decreasing the sample 
Interval. 

A computer algorithm was developed which automates the analysis 
of the large number of sample points encountered In applying this 
method. A listing of the program Is shown In Figure 21. The Input/ 
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e(t) 

5V/div 



t 'lOtis/div 


Waveform excitation into a 1/2 mil 
Supermalloy cut core 


SCL ■ 2. 557x10*^5 ,jp/^,2. 27 ^ 



Internal 

e(t} *V 

dgi/dt ' maxwell/sec 

SCLwattfIb 



1 

0 

0 ^ 

0 



2 

2.5 

25x10® 

.013 

-ive 

half 


3 

6.5 

65x10® 

.110 

A 

4 

9.5 

95x10® 

.261 

cycle 


5 

11.5 

115x10® 

.403 


6 

11.0 

110x10® 

.364 



7 

8.5 

85x10® 

.203 



, 8 

5 

5x10® 

0.000 



1 

0 

0 , 

0 



2 

3 

30x10® 

.027 



3 

6.5 

65x10® 

.157 



4 

8.2 

82x10® 

.266 

+lve 

half 


5 

8.0 

80x10® 

.297 


6 

7.0 

70x10® 

.186 

cycle 


7 

8 

6.5 

5.8 

65x10® 

58x10® 

.157 

.121 



9 

5.0 

50x10® 

35x10® 

.087 



10 

3.5 

.039 



11 

3.0 

30x10® 

.027 



12 

1.0 

10x10® 

.002 



ll3 

0 

0 

0 


Using equation 3.16e 


1.354 +__1J4L_ 

SCLt • 2x8x12.5x10^ 2x13x8. 77id0^ , q watt'pound 

96.9x10“^ 

Figure 20. * An example to demonstrate the application of method 3 
and the specific core loss equations. 
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TJ-0 
TK«K 

FOfl« 1 . / ( 2 . «T I HC*T 1 . E -€ ) 

FOF F« 1 . / ( 2 . *T I NC*TK9 1 . E -« ) 

READ IN ♦ PRIMARY mtTACE VALUES 
X 10 l-l,J 
REAO(S,100)V(1,I) 

OPOT(1,I)-V(l,n/(T«1.E«OI 

CONTINUE 

READ IN > PRIMARY VOLTAGE VALUES 
00 20 l-I^X 
REAO(S,100) V(2,l) 

DPDT(2, I)-V(2,I)/CT«1.0F-8) 

CONTINUE 

OUTPUT OF INITIAL DATA ENTRY 
NPITE(B,195> 

URITE(6,6nO) 

WRITE (6,200) T 
OUTPUT ♦ VOLTAGF INFO 
WRITE(*^,210)FON 
URITE(6,250) 

00 50 l-l,J 

URITF(6,»*0) va, l),DPOT(l, I) 

CONTINUE 

XTPUT - VOLTAGE IfFO 
URITC(6,220)FnFF 
WRirF(6,250) 

X *0 l-l,K 

XITE(6,2*b) V97,I),DXT(2,I> 

COHTIXF 

MATERIAL IOFNTIFICATION ANT TMIPKNESS 
PEAf>(5,115) PMATl 
REAn(5,100lT*Mrx 

if(cmatl.eo.soo»th) X TO 59 

IF(CMAT..cn.$OPEIHf)X TO SX 
|F(CMATL.FO.FERR)GO TO 5*» 
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X TO 60 
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!miTE(6,K30> 
an TO 60 
WRITE(6,520) 

WOITE(6,r.79) 
r.n TO CO 

•l»»TE(6,53n>TH|rx 
tJRITF(G,610) 
r.O TO 6C 

rAiruiAT: rflf’VF^TFp opfratP'O Forr*T?Ty 
FRE0-l./({7r*'C*TJ*TIOC*TE)»l.'’-r.> 
lJRIT6(G,*OniF'>^0 
f>FAP(5,12n>'- 
RFADi 5,1'*9)P 
'TFAP(5,19C)f' 
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P'’ l*3,J 

^rui^ i)-r,#(ppDT(i^ n*.r)*rp9» «*f >) 

5rL9»5ci(i. I)*r.rir 
C9MT.m'F 

Tsci.o-';rLn/(rJ*‘‘.*Fpr) 

:'»'iTF(r,<f>c)T';cio 
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5CL(2, l)-9*fPPP7(-^ l)..r)-fF9cc.*(-* ,) 

^riF-sn (2, 1 ,*S9 lf 
CO riTINtlf 

TSrLE-5CLF/(T9*'».*FnFr) 

’"iTE(r 

C'Lr*.LATr T9T5L SPFriFir ronj- ^9"; 

TSTL- (T59l''*T'-.'‘t c).cren 
’•PITECr, r29)T59L 
For>r,-T{F’»r, 7) 

F0Rr»A7( Z' 

F0'»rA7(Al) 

FPPMATfPK.-^) 

FP^rATMi 1,2* 'l-’pi'T PAta 9'ir ,vArv-|^M 

rOi-”AT{' Ni'r'pr^ or pri'-A'»v t,,/> 

pi'fVTf* PO'^lTIVt V''I.TA9F FP'T'ALfPT orry- ' , r -9 . 7, » M?',//) 

F9i:rv,7(//, ' MFCATP"' V''!TA9' rm*|*'AlF»fT f FOir'-rrv. ' , • M7»,//> 

formak * ♦ I vr v''* TA^F-'m Tr np/’'T-fAv/':'’9*, /) 

rnPMAT(» -1'/- *.,LTA'',F-v9r T*: pp/nT-»*Ay/';pr»,/) 
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ro9**AT(ir»i,/, M7'spp''»^i9 99^r 1,95c r^ri > 59. oRTHOPAt matfrial*) 

F99:<at( u-i, /, 19 ,r rrpF i P5*: (',90 F9r»,cff^*^' hm. :p. pr^MAriPy R9 r-ATCoi.L*) 

Fr)Pr'AT(iri,/,2? *sr r xrf fo" ffpoitf «atf9|ai») 

FOnpAiriHl,/,: ,X'. • 999F (59L) '*''9 * , F ■» . I , ' »*M 'iPPCWALlO'' TC9 1 Al * ' 

coprAKiM’,* f ^ RFRATiHr FrFoi'FNrv* » , C7 9, 7 ^ » h?M 

fo9«at(/,' s^i wATr-r:cr/L9-9Y''L<'' ) 

FO^MATr/,* 59L rPF- ,*'?9.;,* ’'ATT-^fc/' n-9Vrt fM 
Fnp»AT(/,* TOTAL 5f L»’ , f 29, 7, ' /1 9 ‘ , /// ) 

FOomaT(2*:(, 

F0RPAT(rx, 

Fnpp<AT(2lr, •*•*••••#••••*•*) 

’ ) 

FOPMATCPX, ♦•**♦*••••’ ) 

STOP 

FNO 




Rguro 21. - Program listing Ibr SCI analysis program. 
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output formats are shovm in Figure 22. 

The same example shovm in Figu;.e 20 was recalculated using the 
computer algorithm with the sample interval decreased to 2 us. The 
calculated SCL for this example was calculated as 0.142 watt /pound. 
This represents an error of only 1.4 percent when con^>ared with the 
experimentally measured SCL of 0.144 watt/pound. 

This example demonstrates that the third method developed in Sec- 
tion 3.4 may be used to determine the SCL for hi^ly distorted wave- 
forms. Also, by decreasing the sample interval the accuracy of the 
analysis will be improved. 

4 . 3 Examples 

In this section five examples are presented which demonstrate 
the accuracy of the incremental d^/dt method. Three magnetic mate- 
rials were used in the analysis with experimental measurement of their 
SCL mad^ under both sinusoidal and nonsinusoidal excitation. Hie 
waveforms us^d in these analyses are shown in Figure 23. Both ferrite 
3B7 material and Square Permalloy 80 were analyzed using these wave- 
forms. The results of these SCL analyses are shown below: 


Material 

Sinusoidal 

excitation 

Nonsinusoidal excitation 


Calculated 

Measured 

Calculated 

Measured 

Ferrite 

0.388 W/lb 

0.369 W/lb 

0.244 W/lb 

0.247 W/lb 

Square Permalloy 80 

.777 W/lb 

.716 W/lb 

.889 W/lb 

.831 W/lb 


Tne percentages of error shown for these examples are consistent 
with thos . measured for the previous examples. Therefore this method 
may be used for analysis over the wide range of excitation conditions 
encountered in the application of the magnetic materials. 

















OUTPUT POKMAT 

INPUT DATA FOR MACMETIC WAVEFORM ANALY5IC 


NUMBER OF PRIMARY ‘HJRNS- 0, lOOOOOOE 02 


POSITIVE VOLTAGE EQU I VALENT FREQUENCY- 


O.IOOOOOIE 04 MZ 


♦IVE VOLTAGE-VOLTS OP/OT-MAX/SEC 


IMPUT PORMAl 


10 
12 
12 

2. , 

. ? 6 ! 
.765 
1.218 
1.587 
1.848 
1.983 
1.983 
1.848 
1.587 
1.218 
.765 

.261 J 
.261 " 
.765 
1.218 
1.587 
1.848 
1.983 
1.985 
1.848 
1.587 
1.218 
.765 
.261 . 
s 

1 . 


0 of turns 
0 of 4-lve Maple PTS 
0 oi *lve eeaple PTS 
At In veee 


We(t) veluee 


y.«(t) veluee 


ae ter lei type 
tepe thickneee 


4.430^-13 s “ eonetent 

2 . 

.825 


:} 


eap<meati 


0. 1763999E 00 
0. 5132999E 00 
0. 7995999E 00 
0.1007000E 01 
0.1117000E 01 
0.1117000E 01 
0. 1007000E 01 
0. 7995999E 00 
0.S132999E 00 
0.17C3999E 00 


0.17C4000E 07 
0.5133001E 07 
0. 7996002E 07 
0.1007000E 08 
0.1117000E 08 
0.1117000E 08 
0.1007000E 08 
0.T006002E 07 
0.S133001E 07 
C.17G4CC0E 07 


NEGATIYi VOLTAC; EQUIVALENT FREQUENCY- 


O.IOOOOOIE 04 HZ 


IVE VOLTAGE-VOLTS 

0. 17G3999E 00 
n.sn.^ooou ro 
0. 7995999E 00 
n.lij07000E 01 
0.1117000E 01 
0.1117000E 01 
0.1007000E 01 
0.7995999E 00 
0.5132999E 00 
0.1763999E 00 


SPECIFIC CORE 


DP/DT-MAX/SEC 

0. 17b409C: . 7 

G.S133001E 07 
0.799G002E 07 

0.inO7rA(.£ 08 

0.1i:7000E 08 
0.1117000E 08 
0.1C070G0E 01 
0. 7996002E 07 
0.5153001E 07 
0. 1764000E 07 


LOSS (SCL) FOR O.SE 


00 MIL SUPERMALLOY MATERIAL 


CONVERTER OPEiMTIi^G FREQUENCY* O.IOOOOOIE 04 HZ 

SCL ON - 0.C189252E-0S UATT-SEC/LB-CYCLE 

SCL OFF- 0.6189252E-05 WATT-SEC/LB-CVCLE 


ON 

N> 


TOTAL SCL- 0.1237851E-01 WATT/LB 


Figure 22. - lnput*output formats for automated calculation of SCL of a jfOgneUc material. 
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This method may also be used with the square wave SCL character- 
istics to determine the SCL characteristics for sinusoidal excitation. 
Figure 24 shows the results of an analysis of a 1000 Hz sine wave for 
various levels. The calculated results shown as a bold line are 

compared with the measured sinusoidal characteristics for 1/2 mil 
Supermalloy material. The results show that an extremely accurate 
SCL characteristic may be achieved by this method of analysis. 



Specific core ioss, watt/pound 
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Figure 24 - Comparison of calculated and measured results for specific core loss of 
1/2 mil cut Supermalloy under sinusoidal excitation at 1000 Hz. 



CHAPTER V 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE RESEARCH 
5 . 1 Conclusions 

The research covered in this dissertation provides the power elec- 
tronics designer with a unique analytical tool for investigation of the 
core loss characteristics of magnetic materials* Analysis methods pre- 
viously available were based strictly on sinusoidal excitation volt- 
ages, Core loss analysis has been based on the assumption that eddy 
current and hysteresis loss components are separable. By using a 
novel excitation technique employing a pulse width modulated square 
wave, a true separation of these two loss components was shown experi- 
mentally, This approach provides a straightforward analysis technique 
to characterize the basi : core loss process of new materials (Ref, 27), 
For fixed conditions the width of the B-H loop was shown 

to be a direct function of d4>/dt. If this concept is expanded, a 
family of B-H curves for various B^^ levels may be constructed. 

Using this family, an operating B-H locus may be constructed for any 
type of waveform excitation. This approach was used to explain the ob- 
served difference in B-H width when excited from either a voltage 
source (narrow B-H width) or a current source (wider B-H width). 

The flux characteristics under these two source conditions directly 
influence the observed B-H shape. For sinusoidal voltage excitation 
the flux is cosinusoidal in shape and has variable d(^/dt values over 
one cycle. For sinusoidal current excitation the flux approaches a 
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square wave. This excitation causes correspondingly higher d^/dt 
values than those encountered for cosinusoidal flux conditions. 

These results demonstrate the advantage of using square wave ex- 
citation as an analysis tool to determine core loss characteristics. 

In addition, with square wave excitation the application of Faraday's 
law becomes straightforward as d4>/dt is a constant value. The spe- 
cific core loss characteristics were experimentally measured by hold- 
ing d(|)/dt constant and varying frequencies to develop a family of 
core loss characteristics. 

In Qiapter II the state-of-the-art instrumentati> > system is 
discussed. This sysfem was used to obtain large quantities of core 
loss data quickly and accurately. Through the use of a high power 
operatiot amplifier, and three different feedback conditions, cur- 
rent, voltage, and flux sources were easily simulated. This allowed 
the magnetic characteristics of a large number oi source and excita- 
tion conditions to be compared and measured easilv. With square wave 
excitation d<p/dt rates could be calculated d’ i from excitation 

voltage levels. A true rms voltmeter was used tr measure the excita- 
tion voltage level. Accurate core loss was measured for a large num- 
ber of data points using a wideband wattmeter. This wattmeter has 
the capability of measuring the power in a 100 kHz square wave to 
within 2 percent. Instantaneous energy measurements were obtained 
with a multiplying oscilloscope. A four channel oscilloscope with 
both time base and X-Y capabilities allowed simultaneous observation of 
e(t), i(t), p(t), 4>(t) , and the B-H loop char ac '.eristics of the par- 
ticular magnetic material under test. 
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Also in this chapter four different ruethods for determining stack- 
ing factor were covered. Stacking factor error in the published liter- 
ature is as great as 40 percent. With the measurement methods used in 
this chapter to determine stacking factor, results consistent to within 
5 percent were obtained. Of the four methods used, the oil displace- 
ment method yielded the most conservative value for the stacking factor 
parameter. 

Using the SCL characteristics determined under square wave exci- 
tation, and the concept of separation of loss variables at constant 
frequency presented in Chapter III, three analysis techniques were 
developed which predict the specific core loss of a magnetic material 
under square wave excitation. The three methods calculated this core 
loss on the basis of maximum flux density in the core, rate of change 
of flux within the core, or upon the incremental rate of change of flux 
in the core. 

In Chapter IV the third core loss method was extended to Include 
the analysis of both sinusoidal and nonsinusoidal excitation condi- 
tions. To implement the method quickly and reliably, generalized spe- 
cific core loss equations of the various magnetic materials were de- 
rived using experimental data. Once these equations are derived, the 
specific core losses for various excitation conditions may be readily 
calculated using a hand calculator. For the examples analyzed in this 
dissertation a computer algorithm was written to provide a high speed 
analysis tool which implements the concepts of method three. Analyti- 
cal core loss predictions compared well with experimentally measured 
core loss data. These analysis methods are applicable to a broad 
range of excitation conditions. In addition, the basic SCL cbarac- 
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teristics of any desired waveform excitation may be derived. 

5.2 Applications in Power Electronics 

Most converters in use today excite the magnetic components with 
some well defined flux waveshape. The emergence of higher power level 
requirements has led to the development of certain converters using 
L-C resonant power stages. In these converters the magnetic devices 
are excited by pseudo-sinusoidal currents. The analysis methods de- 
veloped in this dissertation are directly applicable to the analysis 
of the magnetic material behavior under these excitation conditions. 

Since all DC-AC/DC converters operate from some steady-state 
source of raw power the magnetic exciting voltage as a function of 
time is a symmetrical or asymmetrical rectangular wave. The exact 
voltage pulse width is determined by the duty-cycle control mechanism 
which is incorporated into the converter. In terms of timing imple- 
mentation , the available duty cycle control schemes are: 

Constant on-time T , variable off time T_ 

n I 

Constant T^, variable T^ 

Constant (T^ + T^) , variable individual T^ and T^ 

Variable T^, T^, and (T^ + T^) 

In the examples shown in Chapter III, the duty cycle directly affects 
the SCL of the magnetic material. Therefore, **he magnetic losses 
associated with these converters are not fixed, but become direct 
functions of duty cycle. The three SCL analysis methods developed 
allow detailed performance calculations of the magnetic loss experi- 
enced in these converter applications. 

5.3 Optimization of Tape Thickness 

Using the SCL characteristics of a specific magnetic uiaterial 
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the core loss dependence as a function of tape thickness was investi- 
gated. Holding and d^/dt constant for the four tape thick- 

nesses measured, it appeared that an optimum tape thickness existed 
for each of die three magnetic materials tested. 

Previous work (Ref. 28) has shown a similar optimization to exist 
with respect to grain size of the magnetic material. Since the mag- 
netic materials examined in this dissertation were picked from random 
material batches with uncontrolled grain size, additional controlled 
eiqieriments us^.ng the square wave SCL characteristics as a measure 
should be performed to examine this observed condition. If a true 
optimum tape thickness for a given materied exists, the traditional 
practice of using thinner magnetic tapes for higher frequencies would 
have to be reevaluated. 



APPENDIX A 


SUMMARY OF SPECIFIC CORE LOSS DATA FOR FERRITE, SUPERMALLOY, 

AND SQUARE PERMALLOY 80 CORE MATERIALS 
The magnetic core loss characteristics are shown in this appendix 
with SCL (watts/lo) as the dependent variable, and frequency (Hz) as 
the dependent variable. Flux density (gauss) and d4i/dt (maxwells /sec) 
are shown as parameters. While the units given for "CL are not in 
international units, 'they are consistent with those used in published 
data for sinusoidal core loss characteristics. For the tape wound 
materials tape thicknesses of 1/2, 1, 2, and 4 mils are presented. 

The ferrite core loss characteristics are for 3B7 material. 
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Figure A. 6, ■ Specific core characteristics for uncut 1/2 mil SupermaHoy. 
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Figure A. 9. - Specific core loss 











APPENDIX 3 
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d(^/dt 

e(t) 
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e^(t) 

AE. 


eq 


eq 


ex 


i,(t) 

k 

L 


p^(t) 


GLOSSARY OF SYMBOLS 

2 

physical cross sectional area of core, cm 

effective cross sectional area of core, cm' 

flux density, gauss 

maximum flux density, gauss 

rate of change of flux, maxwells -sec ^ 

excitation voltage, V 

core voltage, V 

input voltage to core sample, V 
incremental energy, W-sec-lb ^-cycle ^ 
total energy, W-sec-lb ^-cycle ^ 
frequency, Hz 

positive equivalent frequency, Hz 
negative equivalent frequency, Hz 
excitation frequency^ Hz 
magnetizing force, 
input current to core, A 
transformer equation constant 
mean magnetic path, cm 
negative half cycle increments 
number of turns 

positive half cycle increments 
input power to core, W 
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core 

•^oil 

SCL 

SCL^ 

SCL^ 

SF 

\ore 


Fe 


W 


core 


air 




3 

dansity of core, gm/cm 

3 

density of oil, gm/cm 
flux, maxwells 
maximum flux, maxwells 
Instantaneous flux, maxwells 
specific core loss, W-lb ^ 
incremental specific core loss, W-lb ^ 
total specific core loss, W-lb ^ 
stacking factor 
incremental time interval, sec 

3 

total volume of core structure, cm 
volume of iron, cm^ 
weight of core in air, gm 
hysteresis loss, J-cycle**^ 
displacement weight, gm 



LIST OF REFERENCES 


1. Bozorth, R. M. , Ferromagnetism . D. Van Norstrand Co., Inc., 1951. 

2. Glasoe, G. N. and Lebacqz, J. V., P ulse Generators . McGraw-Hill 

Book Co., Inc., 1948. 

3. Smlt, Jan, Magnetic Properties of Materials . McGraw-Hill Book 

Co. , Inc. , 1971. 

4. Legg, V. E. , ‘‘Magnetic Measurements at Low Flux Densities Using the 

Alternating Current Bridge, “ Bell System Technical Journal , 

Vol. 15, No. 1, Jan. 193t, pp. 39-62. 

Chen, D. Y. , “Comparison of High Frequency Magnetic Core Losses 
Under Two Different Driving Conditions. A Sinusoidal and a 
Square-Wave Voltage, “ I.E.E.E. Power Electronics Specialist 
Conference, June 1978. 

6. Magnetic Inc., D esign Manual Featuring Tape Wound Core s. TWC-3^0R. 

7. Hayt, W. H. , Tr. , Engineering Electromagnetic s. McGraw-Hill 

Book Co., Inc., 1967. 

8. Swift, G. W. , “Power Transformer Core Behavior Under Transient 

Condi t ions, “ IEEE Transactions on Power Apparatu s and Systems , 
Sept. /Oct. 1971, pp. 2206-2210. 

9. Warburg, E. , “Magnetic Investigations,** Ann. Ph ys ik Series 3 , 

Vol. 13, 1881, pp. 141-164. 

10. Brailsford, F. and Fogg, R. , '*Anomalous Iren Losses in Cold 

Reduced Grain-Oriented Transformer Steel at Very Low Frequen- 
cies,** lEE Proceedings , Vol. 113, No. 9, Sept .1* o6, pp. 1562-1564. 


87 



88 


11. Swift, W. M. , Shilliag, J. W. , Bhate, S. K. , and Young, F. J. , 

”Eddy Current Losses in a (110) (OO^ Single Crystal of 3% Si-Fe,” 
IEEE Tran, actf^ -IS on Magnetics , Vol. MAG- 10, ril 1974, pp. 810- 
813. 

12. Shilling, I. W. i^a Houze, G. L. , Jr., "Magnetic Properties and 

Domain Struct ur* in Grain-Oriented 3% Si-Fe," IEEE Transactions 
on Magne ^ , ' v. * . MAG- 10, *’ 0 . 2, Jkt 1974, pp. 195-223. 

13. Pry, R. H. and Beer, C. P. , ’^‘alculation of Energy Loss in Mag- 

netic Sheet Materials Using a Domain Model," doyi i§l_ 2 f_A££li^ 
Physics , Vol. iV, No. 3, Mar. 1958, pp. 532-533. 

14. Halliday, D. and Resnick, R. , Physics for Students of Science 

and Engineering . John Wiley & Sons, Inc., 1965. 

15. Hollitscher, H. , "Core Losses in Magnetic Materials at Very High 

Flux Densities When the Flux is Not Sinusoidal, IEEE Transac- 
tions on Magnetics , Vol. MAG-5, No. 3, Sept. 1969, pp. 642-647. 

16. Mitch, J. E. , Lewis, H. A., and Parnell, R. A., "Production Test- 

ing of Tape Core Materials for Magnetic Amplifiers," AIEE Con- 
ference Paper, June 1954. 

17. Storm, H. F. , Magnetic Amplifiers . John Wiley & Sons, Inc., 1959. 

18. Putnam, T. H. , **Economics of Power Transformer Design," ^EEE Trans- 

actions on Power Apparatus and Systems , Vol. 82, No. 12, Dec. 
1963, pp. 1018-1023. 

19. Lee, F. C. , Yu, Y. , and Triner, J. E. , "Power Converter Design 

Optimization," IEEE Transactions on Aero sp ace and Electronic 
Syst- 13, Vol. ArJ-15, No. 3, May 1979, pp. 344-355. 





89 


20. Alexandersotiy E. F. W. , "Magnetic Properties of Iron at Frequen- 

cies up to 200,000 Cycles," AIEE Journal , Vol. XXX, Part III, 
1911, pp. 2433-2454. 

21. Boon, S. C. and Robey, J. A., "Effect of Domain Wall Motion on 

Power Loss in Grain Oriented Silicon- Iron Sheet," IEEE Proceed- 
ings . Vol. 115, No. 10, Oct. 1968, pp. 1535-1540. 

22. Overshott, K. J. , "The Use of Domain Observations in Understanding 

and Inqiroving the Magnetic Properties of Transformer Steels," 
IEEE Transactions on Magnetics , Vol. MAG-12, No. 6, Nov- 1976, 
pp. 840-845. 

23. Lauers, J. D. , Biringer, P. P. , and Hollltscher, H. , "A Simple 

Method of Estimating the Minor Loop Hysteresis Loss in Thin 
Laminations," IEEE Transactions on Magnetics , Vol. MAG-14, No. 5, 
Sept. 1978, pp. 386-388. 

24. Lauers, J. D. and Blringer, P. P. , "Prediction of Core Losses for 

High Flux Densities and Distorted Flux Waveforms," IEEE Traus- 
actions on Magnetics , Vcl. MAG-12, No. 6, Nov. 1976, pp. 1053- 
1055. 

25. Lauers, J. D. , Biringer, P- P. , and Hollitscher, H. , "Estimation 

of Core Losses When the Flux Wavefon. Contains the Fun imental 
Plus a Single Odd Harmonic Component,’* IEEE Transactions o n 
Magnetics . Vol. MAG-13, No. 5, Spet. 1977, pp. 1128-1130. 

26. Ward, W. P. , Introduction to Electrical Engineering . Prentice 

Hall, Inc., 1952. 

27. Chen, D. Y. , "Amorphous Magnetic Alloys for High Frequency Power 

Electronic Applications," IEEE Power Electronics Specialist 
Conference, June 1979. 



90 


28. Overshott, K. J. and Thompson, J. E. , "Magnetic Properties of 

Grain Oriented Silicon Iron. Part 4 - Relationship of Domain 
Spacing, Grain Size, Sheet Thickness, and Power Loss," Pro- 
ceedings of the lEL , Vol. 117, No. 4, Apr. 1970, pp. 865-868. 



BIOGRAPHY 


Naioe: 

Date of Birth: 
Place of Birth: 


Universities attended: 

Dates 

Degree 

Dates of 
degree 

Cleveland State University 

1965-1969 

B.E.E. 

June, 1969 

Ifaiversity of Toledo " 

1969-1974 

M.S.E.E. 

March, 1974 

University of Toledo 

1974-1979 

Ph.D. 

August, 1979 


Major: Electrical Engineering 

Minor: Engineering Science anr) Mathematics 

Publications: 

1. , ”A Digital Regulated Solar Array Power Module NASA TM 

X-2314 , July 1971. 

2. , with R. P. Iwens and Y. Yu, **Time Domain Modeling and 

Stability Analysis of an Integral Pulse Frequency Modulated DC 
to DC Power Converter,** IEEE Pow er Electronics Specialist Con^ 
ference, 1975 Record, pp. 80-90, IEEE Publication 75 CHO 965-4 AES. 

3. , with Y. Yu, M. Bachmann, and F. C. Lee, ’’Formulation of a 

Methodology for Power Circuit Design Optimization,’* IEEE Power 
Electronics Specialist Conference, 1976 Record, pp. 35-44, IEEE 
Publication 76-CHI 084-3 AES. 



91 






92 


4. , with F. C. Lee and Y. Yu, "Modeling of Switching Regu- 

lator Power Stages With and Without Zero- Inductor-Current Swell 
Time," IEEE Power Electronics Specialist Conference > 1976 
Record , pp. 62-72, IEEE Publication 76-CHI 084-3 AES. 

5. , with I. G. Hansen, "Electric Vehicle Power Train Instru- 

mentation - Some Constraints and Considerations," International 
Electric Vehicle Exposition (lEVE) , April 1977. 

6. , with R. P. Iwens and F. C. Lee, "Discrete Time Domain 

Modeling and Analysis of DC-DC Co »verters with Continuous and 
Discontinuous Inductor Current," Second IFAC Symposium on Con- 
trol in Power Electronics and Electrical Drives, October 1977. 

7. , with F. C. Lee, R. P. Iwens, and Y. Yu, "Generalized 

Computer Aided Discrete Time-Domain Modeling and Analysis of 
DC-DC Converters," IEEE Transactions on Industrial Electronics 
and Control Instrumentation, Vol. IECI-26 , No. 2, May 1979, 
pp. 58-69. 

8. , with Y. Yu, and F. C. Y. Lee, "Power Converter Design 

Optimization," IEEE Transactions on Aerospace and Electronic 
Systems , Vol. AES-15, No. 3, May 1979, pp. 344-355. 

Positions Held: 

1. Teaching Assistant, Cleveland State University, 1967-1969. 

2. Aerospace Engineer, NASA Lewis Research Center, 1969-present. 

Membership in Honorary Societies: 

Eta Kappa Nu 


Pi Mu Epsilon 



